Rapid viral metagenomics using SMART-9N amplification and nanopore sequencing [version 2, peer review: 2 approved]

Emerging and re-emerging viruses are a global health concern. Genome sequencing as an approach for monitoring circulating viruses is currently hampered by complex and expensive methods. Untargeted, metagenomic nanopore sequencing can provide genomic information to identify pathogens, prepare for or even prevent outbreaks. SMART (Switching Mechanism at the 5′ end of RNA Template) is a popular approach for RNA-Seq but most current methods rely on oligo-dT priming to target polyadenylated mRNA molecules. We have developed two random primed SMART-Seq approaches, a sequencing agnostic approach ‘SMART-9N' and a version compatible rapid adapters available from Oxford Nanopore Technologies ‘Rapid SMART-9N'. The methods were developed using viral isolates, clinical samples, and compared to a gold-standard amplicon-based method. From a Zika virus isolate the SMART-9N approach recovered 10kb of the 10.8kb RNA genome in a single nanopore read. We also obtained full genome coverage at a high depth coverage using the Rapid SMART-9N, which takes only 10 minutes and costs up to 45% less than other methods. We found the limits of detection of these methods to be 6 focus forming units (FFU)/mL with 99.02% and 87.58% genome coverage for SMART-9N and Rapid SMART-9N respectively. Yellow fever virus plasma samples and SARS-CoV-2 nasopharyngeal samples previously confirmed by RT-qPCR with a broad range of Ct-values were selected for validation. Both methods produced greater genome coverage when compared to the multiplex PCR approach and we obtained the longest single read of this study (18.5 kb) with a SARS-CoV-2 clinical sample, 60% of the virus genome using the Rapid SMART-9N method. This work demonstrates that SMART-9N and Rapid SMART-9N are sensitive, low input, and long-read compatible alternatives for RNA virus detection and genome sequencing and Rapid SMART-9N improves the cost, time, and complexity of laboratory work.

Rapid viral metagenomics using SMART-9N amplification and nanopore sequencing.

Emerging and re-emerging viruses are a global health concern. Genome sequencing as an approach for monitoring circulating viruses is currently hampered by complex and expensive methods. Untargeted, metagenomic nanopore sequencing can provide genomic information to identify pathogens, prepare for or even prevent outbreaks. SMART (Switching Mechanism at the 5' end of RNA Template) is a popular approach for RNA-Seq but most current methods rely on oligo-dT priming to target polyadenylated mRNA molecules. We have developed two random primed SMART-Seq approaches, a sequencing agnostic approach 'SMART-9N' and a version compatible rapid adapters  available from Oxford Nanopore Technologies 'Rapid SMART-9N'. The methods were developed using viral isolates, clinical samples, and compared to a gold-standard amplicon-based method. From a Zika virus isolate the SMART-9N approach recovered 10kb of the 10.8kb RNA genome in a single nanopore read. We also obtained full genome coverage at a high depth coverage using the Rapid SMART-9N, which takes only 10 minutes and costs up to 45% less than other methods. We found the limits of detection of these methods to be 6 focus forming units (FFU)/mL with 99.02% and 87.58% genome coverage for SMART-9N and Rapid SMART-9N respectively. Yellow fever virus plasma samples and SARS-CoV-2 nasopharyngeal samples previously confirmed by RT-qPCR with a broad range of Ct-values were selected for validation. Both methods produced greater genome coverage when compared to the multiplex PCR approach and we obtained the longest single read of this study (18.5 kb) with a SARS-CoV-2 clinical sample, 60% of the virus genome using the Rapid SMART-9N method. This work demonstrates that SMART-9N and Rapid SMART-9N are sensitive, low input, and long-read compatible alternatives for RNA virus detection and genome sequencing and Rapid SMART-9N improves the cost, time, and complexity of laboratory work.

Molecular characterization and sequecing analysis of SARS-CoV-2 genome in Minas Gerais, Brazil.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in Wuhan, China, is the causative agent of the coronavirus disease 2019 (COVID-19). Since its first notification in São Paulo state (SP) on 26th February 2020, more than 22,300,000 cases and 619,000 deaths were reported in Brazil. In early pandemic, SARS-CoV-2 spread locally, however, over time, this virus was disseminated to other regions of the country. Herein, we performed genomic sequencing and phylogenetic analysis of SARS-CoV-2 using 20 clinical samples of COVID-19 confirmed cases from 9 cities of Minas Gerais state (MG), in order to evaluate the molecular properties of circulating viral strains in this locality from March to May 2020. Our analyses demonstrated the circulation of B.1 lineage isolates in the investigated locations and nucleotide substitutions were observed into the genomic regions related to important viral structures. Additionally, sequences generated in this study clustered with isolates from SP, suggesting a dissemination route between these two states. Alternatively, monophyletic groups of sequences from MG and other states or country were observed, indicating independent events of virus introduction. These results reinforce the need of genomic surveillance for understand the ongoing spread of emerging viral pathogens.

Identification and characterization of the anti-SARS-CoV-2 activity of cationic amphiphilic steroidal compounds.

The ongoing COVID-19 pandemic caused a significant loss of human lives and a worldwide decline in quality of life. Treatment of COVID-19 patients is challenging, and specific treatments to reduce COVID-19 aggravation and mortality are still necessary. Here, we describe the discovery of a novel class of epiandrosterone steroidal compounds with cationic amphiphilic properties that present antiviral activity against SARS-CoV-2 in the low micromolar range. Compounds were identified in screening campaigns using a cytopathic effect-based assay in Vero CCL81 cells, followed by hit compound validation and characterization. Compounds LNB167 and LNB169 were selected due to their ability to reduce the levels of infectious viral progeny and viral RNA levels in Vero CCL81, HEK293, and HuH7.5 cell lines. Mechanistic studies in Vero CCL81 cells indicated that LNB167 and LNB169 inhibited the initial phase of viral replication through mechanisms involving modulation of membrane lipids and cholesterol in host cells. Selection of viral variants resistant to steroidal compound treatment revealed single mutations on transmembrane, lipid membrane-interacting Spike and Envelope proteins. Finally, in vivo testing using the hACE2 transgenic mouse model indicated that SARS-CoV-2 infection could not be ameliorated by LNB167 treatment. We conclude that anti-SARS-CoV-2 activities of steroidal compounds LNB167 and LNB169 are likely host-targeted, consistent with the properties of cationic amphiphilic compounds that modulate host cell lipid biology. Although effective in vitro, protective effects were cell-type specific and did not translate to protection in vivo, indicating that subversion of lipid membrane physiology is an important, yet complex mechanism involved in SARS-CoV-2 replication and pathogenesis.

Interacting Epidemics in Amazonian Brazil: Prior Dengue Infection Associated With Increased Coronavirus Disease 2019 (COVID-19) Risk in a Population-Based Cohort Study.

BACKGROUND: Immunity after dengue virus (DENV) infection has been suggested to cross-protect from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and mortality. METHODS: We tested whether serologically proven prior DENV infection diagnosed in September-October 2019, before the coronavirus disease 2019 (COVID-19) pandemic, reduced the risk of SARS-CoV-2 infection and clinically apparent COVID-19 over the next 13 months in a population-based cohort in Amazonian Brazil. Mixed-effects multiple logistic regression analysis was used to identify predictors of infection and disease, adjusting for potential individual and household-level confounders. Virus genomes from 14 local SARS-CoV-2 isolates were obtained using whole-genome sequencing. RESULTS: Anti-DENV immunoglobulin G (IgG) was found in 37.0% of 1285 cohort participants (95% confidence interval [CI]: 34.3% to 39.7%) in 2019, with 10.4 (95% CI: 6.7-15.5) seroconversion events per 100 person-years during the follow-up. In 2020, 35.2% of the participants (95% CI: 32.6% to 37.8%) had anti-SARS-CoV-2 IgG and 57.1% of the 448 SARS-CoV-2 seropositives (95% CI: 52.4% to 61.8%) reported clinical manifestations at the time of infection. Participants aged >60 years were twice more likely to have symptomatic COVID-19 than children under 5 years. Locally circulating SARS-CoV-2 isolates were assigned to the B.1.1.33 lineage. Contrary to the cross-protection hypothesis, prior DENV infection was associated with twice the risk of clinically apparent COVID-19 upon SARS-CoV-2 infection, with P values between .025 and .039 after adjustment for identified confounders. CONCLUSIONS: Higher risk of clinically apparent COVID-19 among individuals with prior dengue has important public health implications for communities sequentially exposed to DENV and SARS-CoV-2 epidemics.

Neutralisation of SARS-CoV-2 lineage P.1 by antibodies elicited through natural SARS-CoV-2 infection or vaccination with an inactivated SARS-CoV-2 vaccine: an immunological study.

BACKGROUND: Mutations accrued by SARS-CoV-2 lineage P.1-first detected in Brazil in early January, 2021-include amino acid changes in the receptor-binding domain of the viral spike protein that also are reported in other variants of concern, including B.1.1.7 and B.1.351. We aimed to investigate whether isolates of wild-type P.1 lineage SARS-CoV-2 can escape from neutralising antibodies generated by a polyclonal immune response. METHODS: We did an immunological study to assess the neutralising effects of antibodies on lineage P.1 and lineage B isolates of SARS-CoV-2, using plasma samples from patients previously infected with or vaccinated against SARS-CoV-2. Two specimens (P.1/28 and P.1/30) containing SARS-CoV-2 lineage P.1 (as confirmed by viral genome sequencing) were obtained from nasopharyngeal and bronchoalveolar lavage samples collected from patients in Manaus, Brazil, and compared against an isolate of SARS-CoV-2 lineage B (SARS.CoV2/SP02.2020) recovered from a patient in Brazil in February, 2020. Isolates were incubated with plasma samples from 21 blood donors who had previously had COVID-19 and from a total of 53 recipients of the chemically inactivated SARS-CoV-2 vaccine CoronaVac: 18 individuals after receipt of a single dose and an additional 20 individuals (38 in total) after receipt of two doses (collected 17-38 days after the most recent dose); and 15 individuals who received two doses during the phase 3 trial of the vaccine (collected 134-230 days after the second dose). Antibody neutralisation of P.1/28, P.1/30, and B isolates by plasma samples were compared in terms of median virus neutralisation titre (VNT50, defined as the reciprocal value of the sample dilution that showed 50% protection against cytopathic effects). FINDINGS: In terms of VNT50, plasma from individuals previously infected with SARS-CoV-2 had an 8·6 times lower neutralising capacity against the P.1 isolates (median VNT50 30 [IQR <20-45] for P.1/28 and 30 [<20-40] for P.1/30) than against the lineage B isolate (260 [160-400]), with a binominal model showing significant reductions in lineage P.1 isolates compared with the lineage B isolate (p≤0·0001). Efficient neutralisation of P.1 isolates was not seen with plasma samples collected from individuals vaccinated with a first dose of CoronaVac 20-23 days earlier (VNT50s below the limit of detection [<20] for most plasma samples), a second dose 17-38 days earlier (median VNT50 24 [IQR <20-25] for P.1/28 and 28 [<20-25] for P.1/30), or a second dose 134-260 days earlier (all VNT50s below limit of detection). Median VNT50s against the lineage B isolate were 20 (IQR 20-30) after a first dose of CoronaVac 20-23 days earlier, 75 (<20-263) after a second dose 17-38 days earlier, and 20 (<20-30) after a second dose 134-260 days earlier. In plasma collected 17-38 days after a second dose of CoronaVac, neutralising capacity against both P.1 isolates was significantly decreased (p=0·0051 for P.1/28 and p=0·0336 for P.1/30) compared with that against the lineage B isolate. All data were corroborated by results obtained through plaque reduction neutralisation tests. INTERPRETATION: SARS-CoV-2 lineage P.1 might escape neutralisation by antibodies generated in response to polyclonal stimulation against previously circulating variants of SARS-CoV-2. Continuous genomic surveillance of SARS-CoV-2 combined with antibody neutralisation assays could help to guide national immunisation programmes. FUNDING: São Paulo Research Foundation, Brazilian Ministry of Science, Technology and Innovation and Funding Authority for Studies, Medical Research Council, National Council for Scientific and Technological Development, National Institutes of Health. TRANSLATION: For the Portuguese translation of the abstract see Supplementary Materials section.

Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2.

With the emergence of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants that may increase transmissibility and/or cause escape from immune responses, there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant, first detected in the United Kingdom, could be serendipitously detected by the Thermo Fisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69-70, would cause a "spike gene target failure" (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern (VOC) that lack spike Δ69-70, such as B.1.351 (also 501Y.V2), detected in South Africa, and P.1 (also 501Y.V3), recently detected in Brazil. We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all 3 variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675-3677 as the primary target and spike Δ69-70 to differentiate, we designed and validated an open-source PCR assay to detect SARS-CoV-2 VOC. Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence and spread of B.1.1.7, B.1.351, and P.1.

Respiratory Viral Shedding in Healthcare Workers Reinfected with SARS-CoV-2, Brazil, 2020.

We documented 4 cases of severe acute respiratory syndrome coronavirus 2 reinfection by non-variant of concern strains among healthcare workers in Campinas, Brazil. We isolated infectious particles from nasopharyngeal secretions during both infection episodes. Improved and continued protection measures are necessary to mitigate the risk for reinfection among healthcare workers.

Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil.

Cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in Manaus, Brazil, resurged in late 2020 despite previously high levels of infection. Genome sequencing of viruses sampled in Manaus between November 2020 and January 2021 revealed the emergence and circulation of a novel SARS-CoV-2 variant of concern. Lineage P.1 acquired 17 mutations, including a trio in the spike protein (K417T, E484K, and N501Y) associated with increased binding to the human ACE2 (angiotensin-converting enzyme 2) receptor. Molecular clock analysis shows that P.1 emergence occurred around mid-November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.7- to 2.4-fold more transmissible and that previous (non-P.1) infection provides 54 to 79% of the protection against infection with P.1 that it provides against non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness.

Local Transmission of SARS-CoV-2 Lineage B.1.1.7, Brazil, December 2020.

In December 2020, research surveillance detected the B.1.1.7 lineage of severe acute respiratory syndrome coronavirus 2 in São Paulo, Brazil. Rapid genomic sequencing and phylogenetic analysis revealed 2 distinct introductions of the lineage. One patient reported no international travel. There may be more infections with this lineage in Brazil than reported.

Evolution and epidemic spread of SARS-CoV-2 in Brazil.

Brazil currently has one of the fastest-growing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemics in the world. Because of limited available data, assessments of the impact of nonpharmaceutical interventions (NPIs) on this virus spread remain challenging. Using a mobility-driven transmission model, we show that NPIs reduced the reproduction number from >3 to 1 to 1.6 in São Paulo and Rio de Janeiro. Sequencing of 427 new genomes and analysis of a geographically representative genomic dataset identified >100 international virus introductions in Brazil. We estimate that most (76%) of the Brazilian strains fell in three clades that were introduced from Europe between 22 February and 11 March 2020. During the early epidemic phase, we found that SARS-CoV-2 spread mostly locally and within state borders. After this period, despite sharp decreases in air travel, we estimated multiple exportations from large urban centers that coincided with a 25% increase in average traveled distances in national flights. This study sheds new light on the epidemic transmission and evolutionary trajectories of SARS-CoV-2 lineages in Brazil and provides evidence that current interventions remain insufficient to keep virus transmission under control in this country.

First case of SARS-COV-2 sequencing in cerebrospinal fluid of a patient with suspected demyelinating disease.

The association between coronaviruses and central nervous system (CNS) demyelinating lesions has been previously shown. However, no case has been described of an association between the novel coronavirus (SARS-COV-2) and CNS demyelinating disease so far. SARS-COV-2 was previously detected in cerebrospinal fluid (CSF) sample of a patient with encephalitis. However, the virus identity was not confirmed by deep sequencing of SARS-COV-2 detected in the CSF. Here, we report a case of a patient with mild respiratory symptoms and neurological manifestations compatible with clinically isolated syndrome. The viral genome of SARS-COV-2 was detected and sequenced in CSF with 99.74-100% similarity between the patient virus and worldwide sequences. This report suggests a possible association of SARS-COV-2 infection with neurological symptoms of demyelinating disease, even in the absence of relevant upper respiratory tract infection signs.

Importation and early local transmission of COVID-19 in Brazil, 2020.

We conducted the genome sequencing and analysis of the first confirmed COVID-19 infections in Brazil. Rapid sequencing coupled with phylogenetic analyses in the context of travel history corroborate multiple independent importations from Italy and local spread during the initial stage of COVID-19 transmission in Brazil.